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Limit  Probabilities  in  a Multi-type 
Critical  Age-Dependent  Branching  Process 


by 

Howard  J.  Weiner 


I . Introduction  and  Assumptions 
Let 

(1-1)  Z_(t)  " the  number  of  cells  of  type  j at  time  t starting  with 
one  new-born  cell  of  type  i at  t = 0 with  1 < i < m in  an  m-type 
critical  age-dependent  branching  process  described  as  follows. 

At  time  t = 0,  one  newly  born  cell  of  type  i starts  the  process,  for 
some  1 < i < m.  The  cell  lives  a random  time  described  by  a continuous 
distribution  function 


(1.2) 


G.(t),  G.(0+)  = 0. 


At  the  end  of  its  life,  cell  i is  replaced  by  new  daughter  cells  of 
byPe  J2  new  cells  of  type  cells  of  type  m with  probability 


P. 


1J2J3  Jm 


Define  the  generating  functions,  for  £ = (s  ,...,s  ),  ± = ), 


(s  1 s 


(1,3)  h (s.,...s  ) sh  (s)  = E P . s^1...sJm  = E P. .s- 

1 1 m i .v  ii,"*j  1 m . ii- 

(jl*"Jm)  1 Jm  1 X 


Each  daughter  cell  proceeds  independently  of  the  state  of  the  system, 
with  each  cell  type  j governed  by  (t)  and-h.^s). 

Assume,  for  1 + e = (l+e,...,l+e)  and  1^  = m-vectors, 


(1.4) 


hi (1  + e)  < ® for  1 < i < m. 


This  insures  that  all  moments  of  h^(is)  evaluated  at  s = 1 may  be 
computed  by  partial  differentiations  under  the  summation  sign. 

Define,  for  1 < i,  j < mt 


(1.5) 


m.  . = 


9h. (s) 
i — 

9s . 

J 


= h. .(1) 
1 1 1 

S = 1 J 


and  assume 


(1.6)  m„  > 0 all  1 < i,  j < m, 
and  let  the  first  moment  mXm  matrix  be 

(1.7)  M = (m.  . ) . 

By  standard  Frobenius  theory  [2],  there  is  a largest  eigenvalue  in 
absolute  value,  denoted  p,  which  is  positive. 

The  basic  assumption  of  criticality  is  that 

(1.8)  p = 1. 

It  follows  again  from  standard  Frobenius  theory  [2]  that  there  are 
strictly  positive  eigenvectors  u > 0,  v > 0 (no  notational  distinction 
between  row  and  column  m-vectors)  such  that 
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(1.9) 


Mu  = u 


vM  ~ v 

m 

u * 1 = 2 u.  = 1 

• >t  i 

i=i 

m 

n • v s £ u.v.  = 1. 
~ — ,,ii 

i=l 


Let  the  quantity  denoted  by  Q(u)  satisfy 


(1.10) 
and  also 
(1.11) 


2 

1 m m m 3 h.  (1) 

0 < Q (u)  = t-  £ £ £ k — V u 0 u v . < 

2 i-1  je=l  r-l  Ssl  Ssr  1 r 1 


a2h.i  (i) 

> 0 1 < 1,  m . 

J k 


Further  assume  that  each  distribution  function  G (t 

l 


(1.12) 


4+5 

t dG^(t)  < 00  for  some  6 > 0. 


Let 


(1.13) 


0 < Lt  = 


t dG± (t)  < » , 


and  denote 


(1.14) 
Let 

(1.15) 
Let 

(1.16) 


Ft(s,t)  = E 


Zil(t>  Zi2<fc>  Zim(t)‘ 


S1  S2 


in 


Z.(t)  = (zil(t),zi2(t),...,z.m(t)) 


SfCO  = (Nil(t),Ni2(t),...,Nim(t)) 


) satisfies 
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denote  the  m-vector  with  entries 


(1.17)  j (t ) = total  progeny  of  type  j born  by  t in  the  above 
critical  m-type  process  starting  with  one  new  cell  of  type  i. 

Denote 

<1.18)  H^s.t)  = E 


Nil(t)  Ni2(t) 


• s 


N.  (t) 

im 

m 


For  k = (k19...,km)  a vector  of  non-negative  integers,  it  is  to  be 
shown  that  as  t -*  00 

(1.19)  P[Z.  (t)  =k] 

t 

for  k 0,  and 


(1-20)  Pt^Ct)  = kj  ~ c > 0, 

where  the  constants  c may  be  found  by  a recursive  method. 


II.  Integral  Equations  and  Approximations 
Denote 


(2.1) 


m 


P.  (t)  = P[Z.(t)  > 0]  = ( £ u.v.ia^u./QCu), 

i=l 


where  ()  (0, ...  ,0)  and  u.  > v means  term  by  stenn  strict  inequality  for 


vectors,  and 


(2.2) 


See  [3] . 


m m m 


Q(u)  ~ 2 Z S 11  hiir(1)ulUrV-f  < “• 
i=l  r=l  Z=l  lZr  1 r x 
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Denote  by  1 - P(t)  the  m-vector 


(2.3) 


LzU&L  = (i-p  (t)si-p  (t),...,i-P  (t)) 


m 


Theorem  1.  Under  the  assumptions  of  section  I,  for  k = (k.^ , . . . ,k  ) an 


m 


m-veetor  of  non-negative  integers  not  all  of  whose  terms  are  zero,  for 

£ -»  oo 


(2.4) 


PfZ^t)  - k]  ~ ~ 
t 


where  c > 0 is  a constant  depending  on  k,  i. 


Proof.  First  consider  the  case  k = e . , that  is,  k = e . = (0,0,...  0^1 

3 j 

where  the  1 is  in  the  jth  place,  all  other  places  have  a zero. 

From  [3] 


(2.5)  F (s,t)  - s (1-G  (t))  + f h (F(s,t-u))dG.(u), 

^ J.  i.  * _ X 1 


where 


(2.6) 


F(s,t)  = (F1(s,t),F2(s,t)...Fm(s,t)) 


Under  the  assumptions  of  the  Theorem, 


(2.7) 


P.  .(t) 


3Fi(£»t) 


9s . 


= Ptz.(t)  =e  ] 
s=0  ^ J 


Applying  (2.7)  to  (2.5)  yields,  for  6..  the  Kroneker  delta, 

ij  ’ 


(2.8)  P.  (t)  = 6. .(1-G. (t))  + 

ij  1J  1 


e*  t m 


2 h,-  «(£(0»t~u))P«<(t-u)dG . (u), 

0 jft=l  1 


all  1 < i,  j < m. 


..00) 
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The  proof  for  the  . . (t)  will  continue  by  a series  of  claims. 

13 

Claim  I.  For  arbitrary  c > 0,  b > 0,  sufficiently  large  T > 0, 
the  functions  Q. . (t)  defined  for  1 < i,  j < m by 


(2.9) 


for  some  0 < e < 1, 
large. 


cu.  t ^ , t > T 


Qij(fc) 


t < T 


satisfy,  for  all  1 < i,  j < m,  and  all 


t sufficiently 


(2.10) 


m 

S v.Q. .(t)  > v (1-G  (t)) 
i_l  J J 


+ 


»t 
J 0 


m m 

2 v 2 h.  .(F(0) ,t-u))Q.  (t-u)dG  (u)  . 
i-1  3 1 


Proof  of  Claim  I.  Writing  the  integral  term  on  the  r.h.s.  of  (2.10)  as 
pc  / 2 pc 


0 

I, 


t/2 


for  t :»2T,  note  that  by  (2.9),  a Taylor  expansion  o 


_ (2~e)  * (3“6) 

and  a Newton  expansion  of  (t--u)  ' until  the  power  t 


f h. j about 

X AJ 

yields,  as 


(2.11) 


F(0,t-u)  = l-P(t-u) , 


™t/2  m m r~ 


(2.12) 


2 v.  2 
°0  i=l  1J6=1 


m 


m 


hii(l)  - S Pr(t-u)hiJtr(l)  + o(  2 Pr(t)) 


r=l 


r=l 


(t-u) 


2-e 


dG^.  (u)  + o(t  ^ ^) 


for  some  6 > 0,  and  hence  (2.10)  equals, 


to  terms  of  order  t 


(3-e) 

9 


6 


(2.13) 


t/2  m 
S v 
0 i=l 


cu 


i ,,  , (2-e)u  , ,1,  x 

___  q + + 0(-)) 


m m m 


dGi(u) 


m 


■t/2  (.fx  ^'WrW^  f Vl»V 


t<3-£)0(u, 


dG  (u)  + o(t  ^ ) 


using  (1.9),  (2.1),  (2.2). 

(2.14)  l-G^(t)  = o (t  ^ ^)  for  1 < r < m 

(2-15)  J*  udG  (u)  = U.  + o(t  ^+^))  for  1 < i < m 

0 1 ~ ~ 


use  of  (2.14),  (2.15)  in  (2.12),  (2.13)  yields  the  result  of  the  Claim  I. 
Claim  II.  If  for  1 < i,  j < m,  (t)  are  differentiable  and 


(2.16) 


R.  . (0)  = 0 
ij 


m pt  m m 

(2.17)  Ev  R.(t)  > S v j S h (l-P(t-u))R,  (t-u)dG.(u) 

J (<)  o i=i  xje=i  xlL  1 


i«l 


for  all  t sufficiently  large,  then 


(2.18) 


m 


S v.R.  . (t)  > 0 

«) 


i=l  J 


for  all  t sufficiently  large. 


Prop f of  Claim  II.  Expand  h^ (l-P(t-u))  in  a Taylor  expansion  about  h^,  ^ ( ^) , 
The  r.h.s.  of  (2.17)  becomes,  after  splitting  the  integral  into  two  parts 

pt /2  pt 

+ J ’ 

0 t/2 

7 


pt/2  m 

r 111 

m -i 

(2.19) 

h ^(1)  - E Pr(t-u)h  ^ (1)  + o(  2 P (t/2)) 

0 i=l 

L r=1 

J 

Rij (t"u)dGi (u) 


+ o(t-(4+6)). 


which  simplifies  by  (1.9),  and  plugging  back  into  (2.17),  where  we  now  let 

m 

(2.20)  R (t)  = Ev.R  (t), 

J i=l  J 

that,  for  all  t sufficiently  large,  assuming  (say)  the  upper  inequality  in 
(2.17), 

(2.21)  R. (t)  > f R. (t-u)dG. (u) . 

J J o J i 

and 

(2.22)  Rj (0)  = 0. 

First  note  that  from  (2.21),  it  cannot  be  the  case  that  R^(t)  = 0 
on  any  interval  [0,T] . For  if  so,  set  t = T in  (2.21)  to  obtain  0 > 0, 
a contradiction. 

Suppose  now  that  there  is  a first  t >0  such  that  R. (t  ) = 0, 
rr  o jo* 

Rj  (t)  £ 0,  t < tQ.  Set  t = tQ  in  (2.21)  to  obtain  a contradiction  if, 

in  fact,  Rj  (t)  > 0 for  t < tQ.  If  R^.  (t)  < 0 for  t < tQ,  by  continuity 

there  is  a T < t such  that  S R . (T) ] > Ir. (t)  I all  t < t . Set  t = T in 
— o j j o 

(2.21)  to  obtain  an  immediate  contradiction. 

This  contradiction  yields  that  R^.  (t)  ^0,  t > 0.  Hence  if  the  upper 
inequality  of  (2.18)  fails,  it  must  be  that  R^  (t)  < 0,  all  t > 0.  Again, 
for  any  interval  [0,t^],  there  is  a t^  < t^  such  that  ]R (t^ ) | > |R(t) | , t < t^. 
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bet  t — in  (2.21)  to  obtain  a contradiction.  This  suffices  to  prove 
Claim  II,  as  the  argument  for  the  lower  inequality  is  similar. 

Claim  III.  The  functions  defined  for  arbitrary  a > 0,  b > 0,  T > 0, 


(2.23) 


satisfy,  for  1 < i,  j < m,  the  system 


t < T 


t > T 


(2.24) 


R.  .(t)  = 

i j 


f . . (t)  + 


nt  m 

£ h.  „(l-P(t-u))R„. (t-u)dG. (u) 
0 £=1  xXl  1 


where , as  t -»  ra5 


(2.25)  . f ..(t)  = o(t‘3). 

ij 

Proof  of  Claim  III.  The  integral  on  the  r.h.s.  of  (2.24)  with  (2.23) 
substituted  becomes,  after  a Taylor  expansion  about  h.  .(1),  less  than 
or  equal  to 


(2.26) 


*b/2  mm  m bu. 

Z (h  _(!)  - Z P (t-u)h  (1)  + o(  Z P (t/2))) dG.(u) 

0 X=1  r=l  r lXr  r=l  r (t-u)2  1 


+ G(t)  - G(t/2 ) 


and  similarly,  the  integral  on  the  r.h.s.  of  (2.24)  is  greater  than  or  equal 
to  the  integral  term  of  (2.26). 
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The  integral  in  (2*26)  equals 


tn 


(2.27)  b 

tJ 


.t/2 


u 


|(1  +~+o(i)) 


» m hite(i)(  W.VVl 

£ £ ^-3 

i=l  r=l  Q (u)  t 


X dG  (u) 
1 


+ o(t-3) . 


Since 


(2.28) 


f 


■t/2 


udG^(u)  = pL  + o(t  ), 


and 


(2.29) 


G(t)  - G( t/2)  = o(t"4) 


this  suffices  for  the  proof  of  Claim  III. 

The  proof  of  Theorem  I may  now  be  completed. 
Define  the  iterative  sequence  for  n = 0,1,2,.. 


(2.30)  P 


r>t  m 


'0  JN1 


with 

bu. 

(2.31)  P^(t)  = P(0)ij<‘t^  ~Rij^  = ~T’  for  t su£ficiently  large. 


Then,  suppressing  the  i,j  part  of  the  subscript, 


(2.32)  |P(1)(t)  -P(Q)(t)|  = i^i1(l-Gi(t))  - f±i(t)  ! = o(t~3) 


ij 


which  implies  that 
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(2.33) 


for  some  c > 0. 


eu. 


P(l)ij(t)  ~ 2 


Assume,  by  an  induction  hypothesis,  that,  for  some  constant  c (the 
value  of  c may  change  from  expression  to  expression) 


(2.34) 


cu, 


P(flij<t)  ~ ~2 


for  0 < £ < n-1 , 


Then,  by  (1.9). 


rii/2  m 


(2-35)  !P(n)ij(t>-P(n-l)ij<t>!^U0  (t-u)dGi  (u)  ^n_1}ij  (t)  | 


+ 2b(l-G(t/2)). 


Again  by  (1.9)  and  a Taylor  expansion  about  h^(l). 


(2.36)  |P(n)ij<t)-P(n.1)I;|(t>|  < \f_ 


/2  m 

E 

0 &=1 


m 


hi£(i>- 

r=l 


m 

+ o(  E P (t/2 ) ) 
r=l 


P(n-l)£j  (t-u)dG.(u)  - P(n_1);Lj  (t)  | +2b  (l-G(t/2)  ) 


or 


(2. 3?)  ip(n)ij<t)  - P(n-1)ij 


nt/2  m 

(t)i  < |c  S 

o je=i 


m 

hU(i>  - E,rr(t-n)hUr(i) 

t-  r'1 


m 


+ o(  E P (t/2)) 

r=l  J (t-u) 


U£  CUi  . 

dG.(u) f | + 2b(l-G(t/2))  , 


2 i 


where  the  two  constants  labeled  c in  (2.37)  have  the  same  value. 


11 


Hence 


j ^ 

<2-38)  VlHj(t)i  < |cJo  ‘'^a+^+o(i))  + |0(i)|]dG1(u) 


cu. 


+ o(t"3) 


or 


(2.39) 


'P(n)ij (t)  “ P(n-l)ij 


(t)i  < 0(4r) , 


which  implies  that,  for  t -*  00 


(2.40) 
and 

(2.41) 


cu. 

P(n)ij  (t)  ~ * 


cu. 


P(n)ij (t)  < y"  for  appropriate  choice  of  c. 


Claims  I,  II  yield  that  for  arbitrary  c > 0,  one  may  write 


(2.42) 


cu. 


for  all  t. 
Hence 


pt  ra 

(2.43)  |Plj(t)  -P(n)iJ(t)|  < Jq  ^ihi£(l^tzui)|PiJ(t-u)  -P(n.1)jtJ(t-u)ldG. 


(u). 
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Denote 


(2.44) 

%)(t)  S 

Then,  using 

(2 .40) , (2 .41) , (2 .42)  in  (2.43)  yield 

(2.45) 

4<n)<t)  ^ V<n'l)<t"U)‘iGi<U) 

Iterating  (2 

.45),  and  denoting 

(2.46) 

A*G  = A(t-u)dG(u), 

0 

(2.47) 

A(n)(t)  £ VGi(n)(t)  < G-R)(t) 

where 

( \ 

(2„48)  G^  7 (t)  is  the  n-th  convolution  of  G. (t) . 

Let  {X£}  be  a sequence  of  I.I.D.  r.v.s.  each  with  distribution 
function  G^t),  and  let 

(2.49)  S = E X.  . 

n i=l  1 

Then 

(2.50)  G.(n)(t)  = P[Sn  < t]  = P[Sn-nM.  < t-n^.j. 
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By  Chebyshev's  inequality. 


(2.51) 


, v Var S 

Gi  s 


Kn 


(n|i.  - t)  (n  |i.  - t)' 

i 1 


2+e 


Let  n > t . Then  (2.47)  yields,  for  each  t, 


(2.52) 


■(2+e) 


Then  (2 . 52 ) , (2 .40) , yield  that 


(2.53) 


pij(t)  ~7 


for  some  c > 0. 


Hence  the  theorem  is  proved  for  k = e . 

In  general,  note  that  for  k = (k.  , . . . ,k  ) and  T,  k. 

1 ” i=l  1 

the  degree  of  differentiation. 


= k,  where  k is 


(2.54)  Pk(t)  ^ 


i ak  i 

7,k  ,...kV  “X =r[2.(t)-y. 

11  m os-  • • • Os  ~ ~ 

I m 


Then,  successively  differentiating  (2.5)  according  to  (2.54)  will 
yield,  by  Leibniz'  rule  for  successive  differentiation 


(2.55)  Plk(t)  - f (t)  + f ” h0(l-P(t-u))Pffli(t-u)dG1(u), 

0 X—  1 — 


where  f^(t)  consists  of  sums  of  terms  of  the  form 
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(2.56) 


/A  m 


0 i,r, 


...( 


— =r h . (F (s , t-u) ) 

os .ds  •••  1 — v— ’ " 

x r 


s=0 


P%(t'u)Prb(t"u)“'dGi(u) 


where 


(2.57)  p < k and  the  sum  of  degrees  of  q,  b,‘“  total  less  than  k. 


By  induction,  assume  that  for  degree  of  k < n 


(2.58) 


Pik(t) 2 for  1 £ 1 < m 


where  c > 0 may  change  in  value. 

By  (2 . 55) - (2 . 58) , (2.55)  may  be  written 


(2.59)  P.k(t)  - f.k(t)  + 


f±  m 


where 


(2.60) 


fik(t)  = °<fc  > 


Hence  the  techniques  of  section  2 used  in  the  proof  of  the  theorem 
for  P^(t)  may  be  used  again,  since  (2 .59) , (2 .60)  is  of  the  form  (2.8). 
Hence  Theorem  1 is  proved. 

Theorem  2.  Under  the  assumptions  of  section  I,  for  k = (k. , . . . ,k  ) a 

1 m 

vector  of  strictly  positive  integers,  as  t 

(2 -61)  P[N.,(t)  = k]  ~ c > 0 

the  constant  c will  depend  on  k,i  and  can  be  obtained  by  a recursive 

argument . 
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Proof.  Denote  for  1 < 1 < ®f 


(2.62)  H.(s,t)  sh.(s.,...,s  ,t)  = E(s. 

i i I m jl 


"il (t) 


N.  (t) 
xm  . 
s ) . 

m 


The  law  of  total  probability  yields 


(2.63)  ^(s.t)  = s^l-G.  (t)  + 


h. (H(s,t-u))dG.  (u)) 
0 1 “ ~ 


where 


(2.64)  H(j3,t)  = (^(s.t)  ,0^  (s,t)  , . . . sHm(_s,t))  , 


Since 


(2.65) 


3H  (s,t) 


i = P[N. (t)  = e ,] 
j 's=0  3 


and  note  that 


(2.66) 


H.(0,t)  = 0, 


it  follows  immediately  that 


(2.67) 


Q,  , (t)  s 5 (1-G  (t)  +h  (0)G.  (t))  ~ 6.  ,h.  (0) 

xj  1J  1 11  1]  1 


as  t -*  °°. 
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By  an  induction,  assume,  for  vector  k = (k  , ...,k  ) of  degree 
m i m - 

k - S k , for  k < n 
i=l  x 


(2.68)  Qifc(t)  H Pl^Ct)  = k]  ~ c. 


Then  for  degree  r = n+l,  note  that,  by  Leibniz'  rule,  Q (t)  is  the 
sum  of  terms  of  the  form 


(2.69) 


m 


s\.(H(s,t-u)) 


n o ..  „ Ss  .3s  3s  • • ♦ 

0 -&,r,q,  • • *=1  1 r q 


Q (t-U)Q-r  Ti(t~u)Q_o  (t-u)dG.  (u) 


r>  ^ a 
s=0  — 


vqv 


possibly  premultiplied  by  a6  , the  Kronecker  delta,  each  term  Q with 
degree  k < n,  (in  fact  the  sum  of  all  Q-degrees  is  < n)  the  induction 
hypothesis  applies  to  the  Q's  in  the  integrand  of  each  term  (2.69),  as 
H(0,t)  = 0,  it  follows  that 


(2.70) 


9 hi(H(£,t-u)) j 

Ss  .3s  3s  • • • I 
X r q ' s=0 


* hi£rq- • • ^ » a constant. 


This  discussion  and  (2 .69) , (2 . 70)  suffices  to  prove  Theorem  2. 


III.  Remarks  and  Extensions. 

The  analog  of  Claim  III  in  one  dimension  could  have  been  used  to 
simplify  the  argument  of  [4] , and  in  particular,  to  not  require  the 
"logt/t2"  condition  of  [1], 

u. 

The  argument  using  Claims  I,  II  to  deduce  that  P..(t)  < --1  could  have 

t2-e 

been  omitted,  and  it  would  have  sufficed  to  note  that 

^ik^'  — ^t)  H P[Z^(t)  >0]  ~ ^ > hut  Claims  I,  II  make  the  method  self- 

contained  and  show  how  finer  estimates  can  be  easily  obtained. 
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The  arguments  of  this  paper  can  be  easily  extended  to  obtain  asymptotic 


marginal  probabilities  as  follows. 

Theorem  III.  Under  the  assumptions  of  section  I,  for  I a subset  of  the 
integers  {l,2,...,m},  as  t -»  “ 


(3.1)  P[Z.  . (t)  = k. , all  j € I]  ~ P[Z.  (t)  = k] 

J -x  - 


where  k . > 0 for  i € I , 
J 


(3.2) 


- ^kl6lX,k262I’ ' ' * ,km6mI'> 


and 


(3.3) 


6il 


1 if  i € I 


0 if  i £ I. 


Also,  as  t -* 


(3.4)  P[N  (t)  =k^ ; all  j € I]  ~ c > 0. 

Outline  of  Proof.  To  prove  (3.1),  a modification  of  the  Theorem  1 
of  [3]  is  required,  as  given  by  the  following  claim. 

Claim.  For  I a proper  subset  of  the  integers  {1,2, . . . ,m},  as  t -*  °°, 


Outline  of  Proof  of  Claim.  Denote 


(3.6) 


J = 


h 

i'll 


a © • 

(61I,62I’"  ' s6ml') 
-J-I  “ 


Then  (3.5)  yields,  for  £ = 1 

J “i 


<3. 7)  l-P.(t)  - + 


Ja 


a- 


Pn(t-u): 


T (t-u)  )dG.  (u)  . 

m j x x 


The  argument  for  P^t)  in  [3]  pp.  527-529,  may  be  modified  for  p (t) 
by  noting  that  the  basic  idea  of  bounding  P_^(t)  above  and  below  by 
probabilities  of  non-extinction  for  appropriate  one -dimensional  critical 
age-dependent  processes,  with  known  behavior,  may  be  retained.  A sketch 
of  the  modifications  required  is  as  follows.  Let 


(3.8) 


G (t)  = max  G. (t) 
1 < i < m i 


Ga (fc)  = min  G (t). 

1 < i < m 

The  critical  m-type  processes  ^(t),  Z (t)  defined  on  p.  528  of  [3] 
remain  the  same,  both  with  generating  functions  h^  (_s) , 1 < i < m and  with 
progeny  lifetime  distributions  GQ(t),  G^(t)  respectively. 

Then  denoting  by  P.jCt),  P^Ct)  the  quantities  corresponding  to 
Pil(fc)  for  the  Z.(t)  processes  respectively,  it  follows  as  before  that 
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(3.9) 


Pil(fc)  <Pil(t)  <Pi3;(t) 

and  also  (3.7)  holds  for  P^t)  , PiT(t)  with  Gi'(t)  replaced  by  G^(t), 

Gp(t)  for  all  1 < i < m,  respectively. 

In  [3]  p.  528,  referring  to  the  numbering  system  there,  inequalities 
(3.14),  (3.15)  hold  as  before,  and  relations  (3.10),  (3.12),  hold  with  the 
term 

v. 

(3.10)  E ~ (1-G  (t)  s c(l-G  (t)) 

i € J-I  a a “ 

added  to  the  respective  r.h.s.,  and  subtracted  from  the  r.h.s.  of  (3.16), 
and  relations  (3.11),  (3.13)  hold  with 


(3.11) 


" it  H c(l-Gp(t)) 


i e j-i 


added  to  the  respective  r.h.s.,  and  subtracted  from  the  r.h.s,  of  (3.17), 
(3.12)  0 < c < 1. 


The  basic  results  of  Lemmas  2,  7,  8 of  [1]  used  in  [3]  p.  530  remain 
true  if  the  term  l-G(t)  is  replaced  by  (1-c) (l~G(t))  in  equation  (2.3)  p.  385, 
(3.2)  p.  389,  (3.3)  p.  389  of  [1].  That  lemma  2 of  [1]  still  holds  is 
clear  since  the  argument  there  depends  solely  on  monotonicity  considerations 
unaffected  by  this  substitution.  The  new  lemma  7 of  [1]  p.  389  should 
now  be  stated  as: 


If 

(3.13) 


:(t)  - a = (l-c)(l-G(t))  + j x(t-y)dG(y), 
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Then 


(3,14)  = ^ + 0(~) 

t H -4/ 

where 


(3,15) 


S-1  = 


tdG(t) , 

Jo 


and  as  t “♦  00 


(3.16) 


t3  (l-G(t) ) -»  0. 


That  this  follows  is  well-known  or  can  be  seen  by  the  methods  of  this 
paper  as  follows. 

One  may  verify  that 


(3.17) 


f (t)  2 ~ t 
U 


satisfies 


(3.18) 


f (t)  - a + r(t)  + J f (t-u)dG(u) 

0 


where 


(3.19) 


r(t)  = o(-j), 


Then,  using  the  iteration  scheme 


(3.20) 


X(n) (t)  = (l-c)d-G(t))  + 


(n-1) 


(t-u)dG(u) 


it  may  be  verified  that  as  t -»  for 


21 


(3.21) 


(t)  = at/ia  , 


(0) 


(3.22) 


t lX(n)(t)-x(n-l)(t:)l  ‘ 0<1) 


for  all  n,  and 


(3.23) 


£ |x(t)  - x^  (t)  j = o(l) 


for  all  n sufficiently  large  with  respect  to  t.  - 

Then  (3 . 17)- (3 .23)  yield  the  new  lemma  7. 

The  lemma  8 uses  lemma  7^  so  remains  true  since  the  new  lemma  7 has 
the  same  conclusion. 

The  other  results  embodied  in  revised  lemmas  of  [1]  remain  true  for 
reasons  of  monotonicity  or  convexity  which  are  unchanged,  or  by  reasoning 
as  just  above. 

Hence  the  results  of  [3]  pp.  530  ff.  remain  true  with  the  indicated 
modifications  and  (3.5)  holds.  But  then  (2.8)  holds  with  the  argument 


(3.24) 


where 


(3.25)  izZj_-£(t) 


and  all  the  results  of  Theorem  1 go  through  as  before  with  modified  (2.8) 
written 

(3.26)  Pij£(t)  = 6..(1-G.(t))  + JtSh.j0(l-PJ_I(t-u))P£jI(t-u)dGi(u). 
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This  suffices  to  indicate  the  essential  modifications  that  yield  (3.1). 
To  indicate  (3.4),  an  example  will  illustrate  the  general  method.  Let 


(3.27) 


Hij  (t)  = P[N  (t)  = 1] 


Note  that 


(3.28) 


3H (s,t) 


3s . 
3 


s=l -e . 


where 


(3.29)  - (1110111)  with  the  zero  in  the  jth  place, 


Apply  (3.28)  to  (2.63)  to  yield,  noting 


(3.30) 


H(k£j.t)  = 0, 


that 


m 


(3.31)  H..(t)  = 6..h.(0)G,.(t)  + (1-6  )J  E h (0)H  (t-u)dG. (u) , 

J 0 h=l  1 


ij  i l 


A consistent  sequenct  {c_}  satisfying,  for  1 < is  j < 


m 


(3.32) 


m 


Cij  - hj11!'0)  + 


is  constructed,  since  note  that 


Then  the  method  of  Theorem  1 may  be  employed  by  defining  the  iterative 


s cheme 


<3'34>  Hij(n+l>(t>  ‘ Vi®6!111  + j1hU<0)H«j(n)<t-">dGi<u> 


1J  1 1 


with 


(3.35) 


It  remains  to  show  that,  for  t -*  00 


(3.36) 


H.  •/  ,n(0  -H.  . (t) 

1 ij (n+1)  ij (n) 


= o(l) 


for  all  n,  and 


(3.37) 


Ivo-hjd/ol  - °0) 


for  sufficiently  large  n with  respect  to  t.  Thus  (3 . 28) - (3 .37) 
suffice  for  (3.27).  The  general  case  follows  by  a tedious  induction  on 
orders  of  derivatives,  using  Leibniz'  rule  to  deduce  the  general  form  of 
the  integral  equation  for  the  quantity  (3.4)  obtained  by  taking  appropriate 
partial  derivatives  of  (2.63),  following  the  arguments  from  (2.65)  to 
(2.69),  and  will  not  be  given. 

The  method  can  be  applied  to  multi-type  sub-  and  supercritical 
cases  also. 
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